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The catalysis of chiral symmetry breaking in the massless weakly coupled QED in a magnetic 
field is studied. It is shown that the effect is due to the dimensional reduction D ^ _D — 2 in the 
dynamics of fermion pairing in a magnetic field. The dynamical mass of fermions (energy gap in 
the fermion spectrum) is determined. The temperature of the symmetry restoration is estimated 
analytically. 
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Recently a new nonperturbative phenomenon, the catalysis of chiral symmetry breaking in an external magnetic 
' field, was revealed in the massless weakly coupled QED The key observation which led to the new effect in 

On QED came from the much simpler NJL model |^-|^ . As was noticed there the catalysis induced by the magnetic field 
On ' should have a universal (model-independent) character. The roots of the universality lie in the effective dimensional 
', reduction D D — 2 in the fermion pairing which results from the dominance of the lowest Landau level (LLL) in 
O i' the low-energy dynamics. As concrete models, the Nambu-Jona-Lasinio (NJL) model as well as QED in 2-1-1 and 3-t-l 
D ' dimensions were considered p|J^,^-p^. The effect of catalyzing the dynamical chiral symmetry breaking under the 
influence of a magnetic field was extended also to the case of external non-abelian chromomagnetic fields (l^Jl^ and 
' finite temperatures |3, |l3| , p^ , as well as to the supersymmetric NJL model |]l5| . 
T— I Here I shall briefly review the analysis of the effect in the weakly coupled QED4 at zero temperature and then 

consider the question of the chiral symmetry restoration at finite temperature . 
] The Lagrangian density of massless QED4 in a magnetic field is 
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C = —F>^^F^, + -[^,z^>^D^^l^], (1) 
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where the covariant derivative is 



5: D, ^d^-ie{Al-' + A,), (0,--X2,-:ei,0) . (2) 

Besides the Dirac index (n), the fermion field carries an additional flavor index a — 1,2, .. . ,Nf. The Lagrangian 
density (|l|) is invariant under the chiral SUi,{Nf) x SUji{Nf) x C/l+r(1) (here I do not discuss the anomaly connected 
D with the current jsu). As is known, there is no spontaneous chiral symmetry breaking at _B = in the weak coupling 
phase of QED jlj]. However, the magnetic field changes the situation drastically: at B ^ the chiral symmetry is 
broken down to S'C/v(^/) = 5't^L+R(^/) and there appear iV| — 1 gapless Nambu-Goldstone (NG) bosons composed 
from fermion- antifermion pairs. The dynamical mass (energy gap) for a fermion can be defined by considering the 
^ Bethe-Salpeter (BS) equation for NG boson or the Schwinger-Dyson (SD) equation for the dynamical mass 
d ' function . Here I follow the first approach. 

The homogeneous BS equation for the iV| — 1 NG bound states takes the form p7[ |: 

XAsi^^y^P) ^ ^'^ J d'^Xid*yid'^X2d*y2GAAAx^Xl)KAiBi;A2B2{xiyi,X2y2)XA2B2i^2,y2;P)GBiBiy2,y) , (3) 

where the BS wave function = {0\T'ipA{x)'iljB{y)\P', P) , /? = 1, . . . ,Nj — l, and the fermion propagator Gab(x, y) = 
{0\TipA{x)ilJB{y)\0)', the indices A = {na) and B = (mb) include both Dirac {n,m) and flavor (a, 6) indices. 
The BS kernel in leading order in a is : 

KA,Bv,A2B2{xiyi,X2,y2) = -47ria(5Q,a242fci7^i„2 7m2mi^P''(2^2 - X2)S{xi - X2)S{yi - y2) 

+4:nia6a,b,Sb2a2'ln^m,1^2n2'^l^'^i^l " X2)S{xi - yi)5{x2 - y2), (4) 

where the photon propagator is 



*This work is done in collaboration with V.P. Gusynin and V.A. Miransky [1,2,3] 
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I^..(.)~/^^fce-(,,.-A^j^ (5) 

(A is a gauge parameter). The first term on the right-hand side of Eq.(||) corresponds to the ladder approximation. 
The second (annihilation) term does not contribute to the BS equation for NG bosons. Then the equation takes the 
form: 

XAB(a;,y;-P) = -47ra / rf-^XidV^AAi (a;, a;i)(5aiQ27^^„^X^^^^ yi; F)(562bi7,^^„,^5BiB(yi, y)2?,.,.(yi - a;i) , (6) 



where, since the lowest in a (ladder) approximation is used, the full fermion propagator G ab(,x, v) is replaced by the 
propagator S* of a free fermion (with the mass m = Tidyn) in a magnetic field. 

Using the new variables, the center of mass coordinate R = {x + y)/2, and the relative coordinate r = x — y, Eq.(|^) 
can be rewritten as 

Xnm{R,r;P) = -4™ J d'^Rid'^nSnn^ ^R- Rl+ ^ ^ 7rn«2X«2m2 (-Rl, n; f )7m2mi 

S^,„,(^^^~R + Ri^V,,,{~n)exp[~ieir + nrAl''\R~Ri)] x exp[iP(i? - i?i)] . (7) 

Here the function Xnm{R, t\P) is defined from the equation 

X^b(x,2/;F) = {0\Ti;A{x)4,B{v)\P,P) = Af^e'^^^ exp[zer^A-'(i?)] X„™(i?, r; P) (8) 

where A'^ are N'j — 1 flavor matrices (tr(A'^A'') = 25^^\ /3,7 = 1, . . . , — 1). The important fact is that the only 
effect of translation symmetry breaking by the magnetic field is given by the Schwinger phase factor in Eq.(^. Thus, 
Eq.(^ admits the "translation invariant" solution Xnm{,RTf"i P) = xi^'i P)- momentum space, one obtain 

r (fq^(PR^(Pk^(Pkn , ,~(P\\ , Q\\ 

Xnm{p-:P) = -i7ra J — ^ ^ cxp (P^ - qi)R J 5„„ J P|| + ^ , + eA^'^' (R^ ) + ^ J 

X lil,n,Xn,rnAk,Ph';a,rnA^irn(^Pn - ^ , + e A^* (R^ ) - ^^V^.^ -p||,kx -p^~ 2eA-*(Rx)) (9) 

(recall that p\\ = (jp ^p^), px = (p^,p^)). Henceforth I shall consider the equation with the total momentum 0. 

The crucial point for further analysis will be the assumption that mdyn <C ■\/|ei?| and that the region mostly 
responsible for generating the mass is the infrared region with k < ■\/|eP|. As it will be seen, this assumption is 
self-consistent. The assumption allows to replace the fermion propagator Snm by the pole contribution of the LLL: 

5<«.(., . . (-i) d - .v.^sig„(«B,) . (i„, 

and Eq.(^) transforms into the following one: 

p(P||,Px) = ^e"^^Pl / d'A^<fk^Sk\\e-''^'^{l-^^W)r 

X |li^p(fc|| , k J^i±^7''(l - ^^W)V,^{k, - , k^ - A^) , (11) 

where p(p||,p^) = {p\\ — mdyn)x(?')(?'|| ~ Wdyn) and £ = l/^/\eB\ is the magnetic field length. Equation (|l^) implies 
that p{p\\,p±) = exp(— £^p^)iy9(p||), where ip{p\\) satisfies the equation 

y^^^) J ~ "^dyn '^ll ~ ™dyn 
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Here 



p\\) 



d k± exp 



(13) 



Thus, the BS equation has been reduced to a two-dimensional integral equation. 

Henceforth I shall use Euclidean space with = —ik^. Then, because of the symmetry 50(2) x 50(2) x 7-" in a 
magnetic field, one arrives at the following matrix structure for f{p\\)'. 



viP\\) = 75(^ + i7i72-B +^11(7 + ijij2P\\D) 

where A, B, C and D are functions of p|. 

Let me begin the analysis of Eq.(|l^) lay choosing the Feynman gauge. Then 

dx exp(-^22,/2) 



(14) 



-v\\) = i^f^ 



(fc|| 



(15) 



and, substituting the expression ( |14| ) for (fipw) into Eq.([lJ), one finds that B = —A, C = D = 0, and the function 
A(p||) satisfies (henceforth I omit the subscript || in momenta) the equation: 



A{p)^ 



a f d^kA{k) 



27r2 J k^+: 



Myn "'O 



dx e^p{~£'^x/2) 
(k - p)2 + X 



(16) 



This equation has been analyzed in [^||9[]. Particularly, as was shown in (see Appendix C), in the case of weak 
coupling a, the mass function A{p) remains almost constant in the range of momenta < < i/£2 g^^^f^ decays like 
l/p2 outside that region. To get an estimate for mdyn at a 1, the external momentum is set to be zero. Then the 
main contribution of the integral is formed in the infrared region with k^ < l/£'^. The latter validates in its turn the 
substitution A{k) A{0) in the integrand of (|l6|), and finally one comes to the following gap equation 



A{0) ^ -^A{0) 
Ztt 



d^k 



fc2 



Myn "'0 



dxe-Kp{-x£'^/2) 
k"^ + x ' 



(17) 



I.e. 



27r 



dxexp{-ml J^x/2) a , 

log X ~ — log 

a; — 1 47r 



'dyn*^ 



(18) 



The latter leads to the following expression for the dynamical mass: 



^dyn 



exp 



(19) 



where O is a constant of order one. The exponential factor displays the nonperturbative nature of this result. 

In order to avoid any confusion let me mention that the electron self-energy (under the influence of a strong external 
magnetic field) in the standard perturbation theory p^Jl9[ | also reveals the double logarithmic asymptotics. However, 
the chiral limit which of the prime interest here is unreachable in the framework of the perturbation theory. The 
latter explains why the results known for several decades did not lead to what is called here the catalysis of chiral 
symmetry breaking by a magnetic field. 

More accurate analysis which takes into account the momentum dependence of the mass function leads to the result 



mdyn 



OVleSjexp 



(20) 



Notice that the ratio of the powers of this exponent and that in Eq.([T^) is 7r/2\/2 ~ 1.1, thus the approximation used 
above is rather reliable. 

It is worth to note that Eq. ( p^ ) reduces to the two-dimensional (another manifestation of the effective dimensional 
reduction in the dynamics of fermion pairing) Schrodinger equation 
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(-A + m2y„ + y(r))vI/(r)=0, (21) 
where the wave function ^'(r) is defined as follows 



and the potential ^(r) reads 



^(r) = / — 4i^)_e»kr (22) 



Note that this potential has the asymptotics V{y) ~ —2a/nr'^ as r ^ oo. The Schrodinger equation with such 
potentials was studied in |20|. Their result for E{a) agrees with the solution in (pO|). 

As is known, the ladder approximation is not gauge invariant. However, the consideration of the general covariant 
gauge [|| shows that the leading term in ln(m^yj^€^), namely \n{m\^^/'^) ~ —-k^/t^J^ol, is the same in all covariant 
gauges. Thus, in spite of absence of the gauge invariance in the ladder approximation, there are no doubts about the 
existence of the effect. 

Up to now I have considered the ladder approximation in QED in a magnetic field. Let me say several words 
about higher order contributions (such as vacuum polarization). As is shown in ||^, taking into account the vacuum 
polarization results in an equation similar to (^6|) but with a replaced by q;/2. Thus, despite the smallness of a. 
the expansion in a is broken in the infrared region in the model. It is a challenge to define the class of all relevant 
diagrams in QED in a magnetic field. Since the QED coupling constant is weak in the infrared region, this problem, 
though hard, seems not to be hopeless. 

Now let me turn to the case of finite temperature. For introducing the Matsubara technique one have to 
change all the expressions in momentum representation according to the rule: 

k° iojn; uJn = TrT{2n +1). (25) 
Thus, the analogue of the equation (^6|) (with m?'{T) in the denominator instead of rn\^^) reads 

. _ oo oo 

dkA{uOmk) f dxexp{—xi'^/2) 



Ai \ "rr. ST f dkA{ujn,k) f 



+ k2+m^{T)J {uj„ - LU,,,)2 + {k - py 



(26) 



Now if one takes n' = 0,p — in the left hand side of Eq.(26) and puts A{ujn, k) « A(cjo, 0) = const in the integrand, 
he obtains the equation 

a ^ r dk r dxexp{-xe^/2) 

TT ^ J u;l + k^+m^T)J {cj^-cooV + k^ + x- ^ 

n- oo_^ p 

It is easy to check that the gap equation ( p7| ) coincides with the equation (58) in |p^ . 
The sum in (p7|) is easily performed: 



oo oo 



^ a f f dkdxexp[—x£ /2] 



IT 



[(7rr)2 + x- m2(r)]2 + (27rr)2(fc2 + m2(T)) 



\ ^WTmHT) \ 2T j Vk^T^ y 2T J j ^ ' 

In the limit T ^ 0, Eq.(^) reduces to the following one 

oo oo 

a f f dkdx exp[—x£'^ /2] 







Vk^ + x^k^ + (ywT^ + yA?2T^^) ' 



(29) 
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which is just what one obtains from Eq.(p^ after performing the integration over ^4 = —iko- 
The equation for the critical temperature is obtained from ( ^8|) by putting m(Tc) = 0: 







dkdxe^^^^^^ 
[l/4 + a;]2 + /c2 



1/4 + a; 
k 



tanh (tt/c) 



1/4 -a; 



coth 



(30) 



where the change of the variables x — > {2TrTc)'^x and k — > 2iTTck was made. 

By assuming smallness of the critical temperature in comparison with the scale put by the magnetic field, TcH. <C 1, 
it is seen that the double logarithmic in field contribution in Eq.(|30|) comes from the region < x < l/2(7rTc^)^, 
I/tt ^ fc < 00. Simple estimate gives: 



l/2(7rTe^)^ 



a 



dx 



dk 



[1/4- 



1/4 + a; 1/4 -a; 



i/tt 

Then, for the critical temperature, one obtains: 



mdyn{T = 0), 



a 1 2 
— log 

47r ^ 



[2(7rT,^)2 



(31) 



(32) 



where mdyn is given by (|19D. The relationship Tc ~ ?Tidyn between the critical temperature and the zero temperature 
fermion mass was obtained also in NJL model in (2+lj- and (3+l)-dimensions |0,f2|. 



In passing, let me just briefly note that the photon thermal mass, which is of the order of ^/aT |22|, cannot change 
our result for the critical temperature (unlike the vacuum polarization) . As is easy to check, it would lead to the shift 
in X for a constant of the order of a in the integrand of (|l]). However, such a shift is absolutely irrelevant for our 
estimate (^. 

In conclusion, the main result of this paper is the analytic estimate for the dynamical mass and the temperature of 
the chiral symmetry restoration in the weakly interacting QED in a background magnetic field. 
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